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Abstract 

The LHCb Collaboration announced two pentaquark-like structures in the J/'f’P invariant 
mass distribution. We show that the current information on the narrow structure at 4.45 GeV 
is compatible with kinematical effects of the rescattering from XdP to J/'ipp: First, it is 
located exactly at the Xci P threshold. Second, the mass of the four-star well-established 
A(1890) is such that a leading Landau singularity from a triangle diagram can coincidentally 
appear at the Xci p threshold, and third, there is a narrow structure at the Xci p threshold but 
not at the XcoP and Xc2j3 thresholds. In order to check whether that structure corresponds 
to a real exotic resonance, one can measure the process A° —>■ K~XciP- If the Pc(4450) 
structure exists in the Xci P invariant mass distribution as well, then the structure cannot be 
just a kinematical effect but is a real resonance, otherwise, one cannot conclude the Pc(4450) 
to be another exotic hadron. In addition, it is also worthwhile to measure the decay T)!^) —> 
J/ijipp'. a narrow structure at 4.45 GeV but not at the XcoP and Xc 2 P thresholds would 
exclude the possibility of a pure kinematical effect. 
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The observation of many different hadrons half a century ago stimulated the proposal of the 
quark model as a classification scheme [1], and helped to establish quantum chromodynamics 
(QCD) as the fundamental theory of the strong interactions. Since then, hundreds of more hadrons 
were discovered. A renaissance of hadron spectroscopy studies started in 2003, and since then a 
central topic is the identification of the so called exotic hadrons. These are states beyond the 
naive quark model scheme, in which mesons and baryons are composed of a quark-antiquark 
pair and three quarks, respectively. Most of the new interesting structures were observed in the 
mass region of heavy quarkonium, and are called XYZ states (for a list of these particles and 
a review up to 2014, see Ref. [2]). In particular, the A(3872) [3] extremely close to the D^D*^ 
threshold is widely regarded as an exotic meson, and the charged structures with a hidden pair 
of heavy quark and heavy antiquark such as the Zc(4430) [4, 5], Z^(3900) [6, 7], Zf{AQ2Q) [8], 
and (10610,10650) [9] would be explicitly exotic were they really resonances, i.e. poles of the 
S'-matrix. Candidates for explicitly exotic hadrons were extended to the pentaquark sector by the 
new LHCb observations of two structures, denoted as in the J/ipp invariant mass distribution 
with masses (widths) (4380 ± 8 ± 29) MeV ((205 ± 18 ± 86) MeV) and (4449.8 ± 1.7 ± 2.5) MeV 
((39 ± 5 ± 19) MeV), respectively [10]. They were suggested to be hadronic molecules composed 
of an anticharm meson and a charmed baryon [11, 12, 13] the existence of which were already 
predicted in Refs. [14, 15, 16, 17]. They were also discussed as a pentaquark doublet in Ref. [18]. 

Normally, such structures are observed as peaks in invariant mass distributions of certain final 
states, and fitted by using the Breit-Wigner parameterization to extract the masses and widths. 
However, such a procedure is problematic. On the one hand, many of these structures are very 
close to certain thresholds to which they couple strongly. In this case, a use of Breit-Wigner is 
questionable and one needs to account for the thresholds. This can be achieved using the Flatte 
parameterization [19] (a method in this spirit for near-threshold states with coupled channels 
and unitarity was recently proposed in Ref. [20]). On the other hand, not every peak should be 
attributed to the existence of a resonance. In particular, kinematical effects may also show up 
as peaks. Such kinematical effects correspond to singularities of the S'-matrix as well, but they 
are not poles. In general, they are the so-called Landau singularities including branch points 
at thresholds and more complicated ones such as the triangle singularity, also called anomalous 
threshold (detailed discussions of these singularities can be found in the textbooks [21, 22]). The 
observability of the triangle singularity was extensively discussed in 1960s (see Refs. [23, 24, 25] 
and references therein), and recently was used to explain some structures including the ry(1405), 
oi(1420) and (j>{2170) [26, 27, 28, 29, 30, 31]. In fact, there were suggestions that some of the 
Zc and Zi, states were threshold effects [32, 33, 34, 35, 36] and the threshold effects might be 
enhanced by triangle singularities [37]. For a general discussion of 5'-wave threshold effects, see 
also Ref. [38] . Therefore, in order to establish a structure as a resonance, one has to discriminate 
it from such kinematical effects. Indeed, this is possible. As discussed in Ref. [39], a resonance can 
be distinguished from threshold kinematical effects only in the elastic channel which is the channel 
with that threshold. The purpose of this paper is to discuss the possible kinematical effects for 
the narrower structure at 4.45 GeV in the LHCb observations and suggest measurements to check 
whether it is a real exotic resonance or not. 

We first notice that the Pc(4450) structure is exactly located at the threshold of a pair of xa 
and proton, (4448.93 ± 0.07) MeV, and 

Mpx 4050 ) - -Mp = (0.9 ± 3.1) MeV. (1) 

If the angular momentum between the Xci and proton is a P-wave, then the two-body system 
can have quantum numbers = (1/2, 3/2, 5/2)“, compatible with the favored possibilities 
5/2+,5/2“ and 3/2“ [10]. The XciP can rescatter into the observed J/V'P by exchanging soft 
gluons. Two possible diagrams for such a mechanism are shown in Fig. 1, where (a) is a two-point 
loop with a prompt three-body production A° —>■ K~XciP followed by the rescattering process 
XciP J/Ip p^ and in (b) the K~p pair is produced from an intermediate A* state and the proton 
rescatters with the Xci into the J/ip p. We will discuss them subsequently. 

It is worthwhile to notice that the Xci can be produced in the weak decays of the Ab with a 
similar magnitude as that for the J/ip. In the bottom quark decays, the charm quark is produced 
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Figure 1: Two-point and three-point loops for the mechanism of the XciP J/il^P rescattering 
in the decay A° — K~J/ipp . 


via the mediation of the VF-boson. After integrating out the off-shell mediators, one arrives at 
two effective operators for the b —>■ ccs transition: 

01 = [c“7'"(l - 75 )c“][s'^7m(1 - 75 )^^], O 2 = [c“7'"(l - 75 ) 0 ^][s^7m(1 “ 75)^“] , (2) 


where one-loop QCD corrections have been taken into account to form Oi. Here, a,f3 are color 
indices, and they should be set to be the same in O 2 in order to form a color-singlet charmonium 
state. The quark fields, [ 07^(1 — 75)0], will directly generate the charmonium state. A charmonium 
with = 1 like the J/ij) is produced by the vector current, while the axial-vector current 
tends to produce the Xci with = 1++ and the 77 c state with = 0^ . Since the vector 
and axial-vector currents have the same strength in the weak operators, one would expect the 
production rates for the J/ip and Xci are of the same order in b quark decays. Corrections to this 
expectation come from higher-order QCD contributions but are sub-leading [40]. In fact, such an 
expectation is supported by the B meson decay data [2]: 

B{B+ J/ipK+) = (10.27 ± 0.31) x 10-'‘, B{B+ XciK+) = (4.79 ± 0.23) x lO'^ (3) 


Having made these general observations, we return to the discussion of the A° decays measured 
by LHCb. We will first focus on the two-point loop diagram whose singularity is a branch point at 
the Xci P threshold on the real axis of the complex s plane, where and in the following yfs denotes 
the invariant mass of the J/tpp or XciP system. It manifests itself as a cusp at the threshold if the 
Xci P is in an S'-wave. For higher partial waves, the threshold behavior of the amplitude is more 
smooth and a cusp becomes evident in derivatives of the amplitude with respect to s. Since we 
are only interested in the near-threshold region, both of the Xci and the proton are nonrelativistic. 
Thus, the amplitude for Fig. 1 (a) is proportional to the nonrelativistic two-point loop integral 



d^q gVA(g^) 

(27r)3 E — mi — m 2 — q^/{2p,) ’ 


(4) 


where mi 2 denote the masses of the intermediate states in the loop, p is the reduced mass and E 
is the total energy. Here, we consider the case for the P-wave Xci P which has quantum numbers 
compatible with the possibilities of the Pc(4450) reported by the LHCb Collaboration, though 
one should be conservative to take these determinations for granted as none of the singularities 
discussed here was taken into account in the LHCb amplitude analysis. If we take a Gaussian 
form factor, /A((f^) = exp (— 2 (f^/A^), to regularize the loop integral, the analytic expression for 
the loop integral is then given by 


Ga(P) 


" (2^)3/2 \ 4 J 


-2k'^/A^ 
2tt 



(5) 


with k = ij2p[E — mi — m 2 + *e), and the imaginary error function erfi( 0 ) = (2/y7r) e* dt. A 
better regularization method should be applied in the future, but for our present study such an 
approach is fine. 
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Figure 2: Fit to the real and imaginary parts of the Pc(4450) amplitude shown in Fig. 9 in Ref. [10] 
with Eq. (6). The blue curve represents the best fit. It is counterclockwise with increasing the 
Jipp invariant mass from 4.41 GeV to 4.49 GeV, the same range as for the LHCb diagram. 


Using an amplitude with the loop function given in Eq. (5), one can get a peak around the 
XciP threshold. In order to have a more quantitative description of the effect of Fig. 1 (a), we fit 
to the Argand plot for the Pc(4450) amplitude depicted in Fig. 9 (a) in Ref. [10] with an amplitude 

A(a) =iV [& + Ga(E)], (6) 

where 6 is a constant background term which may originate from a direct production of the 
A'“J/'0p, and N is an overall normalization. We fit to both the real and imaginary parts of 
the Pc (4450) amplitude by minimizing the sum of the chi-squared values for both the real and 
imaginary parts. The best fit with a real background term has y^/d.o.f. = 1.75 and is given by 
N = 3144, b = —2.9 x 10“"^ GeV^ and A = 0.16 GeV. With a real background term, the amplitude 
in Eq. ( 6 ) can only be complex when the energy is larger than the Xci P threshold, as is evident 
in Fig. 2. The background is in general complex as a result of the fact that the K, J/ip and p can 
go on shell and many A resonances can contribute to the K p state. One sees from the figure that 
the counterclockwise feature of the LHCb amplitude is reproduced, and the overall agreement is 
good. The absolute value of the amplitude in Eq. ( 6 ) with these determined parameters has a 
narrow peak around the XciP threshold as shown in Fig. 3 (a). We have checked that using a 
different form factor A"^/ -|- A^)^ gives a similar result. In both cases, the peak is asymmetric 

unlike the Breit-Wigner form. 

There can be further enhancement around the Xci P threshold due to the presence of nearby 
triangle singularities, also called leading Landau singularities of a triangle diagram, from Fig. 1 (b). 
The leading Landau singularities for a triangle diagram are solutions of the Landau equation [41] 

1 + 27/12 2/23 2/13 = y?2 + 2/23 + yi3: (7) 

where yij = (m| -|- — p^^) / (2 rrii rrij) with mi{i = 1,2,3) masses of the intermediate particles, 

and pij = Pi + pj being the four momentum of the ij pair. To be specific, we let mi, m 2 and m 3 
correspond to the masses of the A*, J/i/i and proton, respectively. Then, PI 2 = p\^ = 

and P 23 = s is the invariant mass squared of the J/ip P pair. It is easy to solve this equation 
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Figure 3: Absolute values of amplitudes in arbitrary units: (a) is for the amplitude in Eq. (6) 
fitted to the Argand plot; (b) is the for the triangle loop integral with the Xci P vertex in a P-wave. 
In (b), we assume the A(1890) with a mass of 1.89 GeV is exchanged in the triangle diagram. The 
solid, dashed and dotted lines correspond to a width of the A(1890) of 10, 60 and 100 MeV, in 
order. 


for any given variable. We solve it as an equation of s, which has two solutions. For an easy 
visualization, we plot in the left panel of Fig. 4 the motion of the solutions in the complex ^/s 
plane. As discussed in 1960s, see e.g. Ref. [23], only one of the singularities can have an impact 
on the amplitude in the physical region defined on the upper edge of the real axis on the first 
Riemann sheet of the complex s-plane, and it is effective only in a limited region of one of these 
variables. Here we want to investigate in which values the A* mass can take so that there can be 
an evident singularity effect in the J/ipp invariant mass, y^. According to the Coleman-Norton 
theorem [42] , the singularity is in the physical region only when the process can happen classically, 
which means that all the intermediate states are on shell, and the proton emitted from the decay 
of the A* moves along the same direction as the Xci and can catch up with it to rescatter. Let us 
start from a very large mass for the A* so that it cannot go on shell in Fig. 1 (b). Decreasing this 
mass, when it has a value 

™i,high = \f^2 - , (8) 

it can go on shell. At this point, the Xci is at rest in the rest frame of the decaying particle A;,, 
and the proton emitted from the decay A* —> K~p can definitely rescatter with the Xci classically. 
This is the point shown as a filled triangle with My. =2.11 GeV, labelled as A, on the solid curves 
in Fig. 4. If we decrease mi further, the Xci will speed up and the proton will slow down. Thus, 
the lower bound of mi for the rescattering process that happens classically is given by the case 
when the Xci and the proton are at a relative rest, i.e. when the Xci P invariant mass is equal to 
their threshold. Thus, at this point the triangle singularity coincide with the normal threshold, 
and one gets 

TOi.iow = 1/-;-m 2 m 3 . (9) 

y m2 + m3 

If mi is smaller than mi^iow, the proton would not be able to catch up with the Xci and the 
triangle diagram can only be a quantum process. For the case of Fig. 1 (b), mijow is given by 
Mji* = 1.89 GeV, labelled as B and also shown as a filled triangle in Fig. 1. In the left panel of 
Fig. 4, in order to move the singularity trajectories away from the real axis, we give a 10 MeV 
width to the A*. For a vanishing width, the solid and dashed trajectories would pinch the real 
axis at mi = mi^high- We can now know on which Riemann sheet of the complex s-plane the 
singularities are located. Since only when mi is between mijow and mi^high (the part between 
the two filled triangles in the figure), the process can happen classically and the singularity can 
be on the physical boundary (if the A* width vanishes), we conclude that the singularity shown 


5 










ReV^ [GeV] 



Figure 4: Left: Motion of the two triangle singularities in the complex plane of ^/s = = 

Mj/^p with respect to changing the mass of the exchanged A* baryon (several values are labeled 
in the plot in units of GeV). In order to distinguish the trajectories from the real axis, we put 
a small imaginary part, —5 MeV corresponding to a width of 10 MeV, to M\-. Only the part 
between the two filled triangles, labelled as A and B, has a large impact on the physical amplitude. 
The thick solid straight line represents the unitary cut starting from the Xci P threshold. Right: 
The corresponding Dalitz plot which shows the region between A and B. 


as the solid curve is always on the second Riemann sheet. On the contrary, the singularity whose 
trajectory is shown as the dashed curve in the left panel of Fig. 4 is on the second Riemann sheet 
when it is above the real axis, and it moves into the lower half plane of the first Riemann sheet 
otherwise. Thus, it is always far away from the physical boundary, and does not have any visible 
impact on the physical amplitude. For an easy visualization of the kinematical region between A 
and B, we show the corresponding Dalitz plot in the right panel of Fig. 4. 

An intriguing observation for the case of interest is that within the range between 1.89 GeV and 
2.11 GeV, there is a four-star baryon A(1890) with 3/2+. Taking = 1.89 GeV, the triangle 
singularity is just at the XciP threshold which can provide a further threshold enhancement. ^ 
Giving a finite width to the A(1890), the singularity moves away from the real axis into the 
lower half plane of the second Riemann sheet (it is located at (4447 — iO.2) MeV for Ma* = 
(1.89 — i0.03) GeV), and the enhancement is reduced. The A(1890) has a relatively small width 
(60 to 100 MeV [2]) so that there can still be an important enhancement. In Fig. 3 (b), we show 
the absolute value of the triangle loop integral with the Xci p in a P-wave for three different widths 
(for a discussion of the triangle singularities in nonrelativistic triangle loop integral, see Ref. [43]). 
There is clearly an enhancement nearby 4.45 GeV even when the width is taken to be 100 MeV. 

In the above, we have shown that kinematical effects can result in a narrow structure around 
the XciP threshold in the J/'i/'P invariant mass of the A[] —>■ K~J/%l>p decay. Gonsequently, a 
natural question is whether such an effect happens at other thresholds, in particular those related 
to the XciP through heavy quark spin symmetry (HQSS). As a result of the HQSS, the operator 
for annihilating a Xci and creating a J/ij: is contained in 

\ {Jh^) = + 4K , ( 10 ) 

where the fields J = 4■ ^ and x = o'! Xc 2 ~ [44, 45] annihilate 

^The mechanism of enhanced threshold effect due to the triangle singularity was recently discussed for the case 
of Zc and states [37]. 


6 









the S-wave and P-wave charmonium states, respectively. This means that the rescattering inter¬ 
action strength for Xc 2 P J/ip p or XcoP J/ip p is of similar size as that for the XciP J/'P’P- 
One might naively expect enhancements at both the Xc 2 P and XcoP thresholds in the J/ipp in¬ 
variant mass as well. However, this is not the case. As we have shown in Eq. (2), at leading order 
in Us, the charmoium is produced by the [c 7 ^(l — 75 ) 0 ] current. This current has no projection 
onto the Xco or Xc2- The production of the Xc0,c2 in the b decays can come only from higher-order 
QCD corrections which are suppressed. Indeed, there is no enhancement at the Xc 2 P and XcoP 
thresholds in Af, decays, which is consistent with our expectation. 

The above analysis is applicable to any b quark decay in which the Xc0,c2 is directly generated 
by the weak interaction. But it would be different if the initial decay heavy particle contains a 
charm or anticharm quark in addition to the bottom quark. Processes of this type include the 
decays of the Be meson and the doubly-heavy baryon Sf,c. An explicit calculation of Be decays 
[46] indicates that the Xc0,c2 can have similar production rates with the Xci- Considering the large 
amount of data on the Be to be accumulated by the LHCb collaboration [47], it appears very 
promising to investigate the XcoP and Xc 2 P threshold effects in the future. In addition, one can 
study the threshold effects in the prompt production of the J/ipp at the LHC, or in the T(IS') 
decays into the XcJPP and J/ip pp. 

In conclusion, what we have shown here is that the present information on the narrow structure 
around 4.45 GeV observed by the LHCb Collaboration is compatible with kinematical effects 
around the XaP threshold: First, it is located exactly at the XaP threshold. Second, the mass of 
the four-star well-established A(1890) coincidentally makes the triangle singularity on the physical 
boundary located at the Xci P threshold, despite a small shift into the complex plane due to the 
finite width of the A(I890), and third, the Xcij instead of the XcO or Xc 2 j can be easily produced 
in the weak decays of the Af, by the V — A current so that there can be an evident effect at the 
XciP, but not the XcoP or Xc2P, threshold. 

Therefore, the most important question regarding the structure around 4.45 GeV is whether it 
is just a kinematical effect or a real resonance. As discussed in Ref. [39], kinematical singularities, 
including both the normal threshold and the triangle singularity, cannot produce a narrow near¬ 
threshold peak in the elastic channel, which is the Xci P in this case. The reason is the interaction 
strength in the elastic channel controls the threshold behavior, and there can be a narrow near¬ 
threshold peak only when the interaction in the elastic channel is strong enough to produce a 
pole in the S'-matrix which corresponds to a real resonance. On the contrary, one cannot simply 
determine the interaction strength for the inelastic channel (xciP J/ip p in our case) because it 
can always interfere with a direct production of the final state. Thus, the question can be answered 
by analyzing the process A[] —>■ K~XciP' if there is a narrow structure just above threshold in the 
Xci P invariant mass distribution, then the structure cannot be just a kinematical effect and calls 
for the existence of a real pentaquark-like exotic resonance, otherwise, one cannot conclude the 
Pc(4450) to be another exotic hadron. 
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